Table of contents.html¶
- ma1 Jupyter
- Using Jupyter
- Cell manipulation
- Configuration
- IPython
- Completion and help
- Shell under IPython
- Magic commands
- Using Jupyter
- ma1 np01 Numpy Introduction
- NumPy - N-dimensional Array manipulations library
- Creating an array
- By redefining its shape
- Mixing values
- Basic Operations
- Browse an array
- Think vector
- Creating an array
- NumPy - N-dimensional Array manipulations library
- ma1 np02 Filtres
- Filter by indices
- Logical filters
- A filter = a logical condition
where
to handle non-filter values- Update a table with a filter
- ma1 np03 Manipulations
- The axes
- Arranging a table
- Reorder axes
- Changing the order of array elements
- Aggregation
- Concatenation
- Stacking
- Splitting
- From Python to Numpy
- Pandas too
- ma1 np05 Notation Einstein
- Introduction to Einstein Notation
- Practical Application
- ma1 np06 Linalg pour le calcul matriciel
- Linalg (linear algebra)
- Basic operations
- Extraction
- Matrix operations
- Linalg (linear algebra)
- ma1 np90 petits exercices
- Numpy - Exercises
- Square Matrix
- Vector Norm
- Sub-Matrix
- Random Vector
- Trace
- Matrix of Multiples of 3
- Count of 9s
- Column with the Smallest Average
- ChessSum
- 2 Minimums
- Rows in Order
- Unique Values
- Magic Tensor
- Tensor Slices
- Numpy - Exercises
- ma20 Rappels sur les matrices
- Vector
- Matrices and linear maps
- Determinant of a matrix
- Standards
- Norm of a vector
- Norm of a matrix
- Properties
- ma21 Transformations isometriques
- Isometric transformations
- Rotation matrix centered at (0,0)
- Properties
- Axial Symmetry
- Translation
- Exercise 1.1
- Rotation matrix centered at (0,0)
- Isometric transformations
- ma22 Changement de repere
- Matrice de passage
- Vecteurs dans le nouveau repère
- Matrice de passage vue comme une transformation
- Points dans le nouveau repère
- Notre souris dans le nouveau repère
- Exercice -- Et l'inverse ?
- Une application linéaire transposée dans le nouveau repère
- Matrice de passage
- ma24 Vectors propres
- $A^n \textbf{x}$
- Vecteurs propres et valeurs propres
- Le cas des matrices de rotation
- Symétrie axiale horizontale
- Diagonalisation d'une matrice
- $A^n \textbf{x}$
- ma25 Drones -- Exercice
- Drone show
- Figure 1
- Figure 2
- Figure 3
- Drone show
- ma26 Vecteurs propres -- Exercices
- Cas d'utilisation des valeurs et vecteurs propres
- Fibonnacci
- Google page rank
- Approche itérative
- Un autre approche
- Cas d'utilisation des valeurs et vecteurs propres
- ma30 ACP
- Principal component analysis (PCA)
- A cloud of dots
- Covariance matrix
- A cloud of dots
- Principal component analysis (PCA)
- ma31 Système d'équations
- Systèmes matriciels
- Résolution d'un système matriciel
- Méthode du pivot de Gauss
- Complexité du pivot de Gauss
- Décomposition LU (Lower, Upper)
- Gauss Jordan
- Comparaison de la vitesse de méthodes
- Erreurs d'arrondi
- Solution au problème d'arrondi dans le cas du pivot de Gauss
- ma32 Conditionnement d'une matrice
- Conditionnement d'une matrice
- Pourquoi ?
- Perturbons la matrice
- Propriétés
- Préconditionnement
- Conditionnement d'une matrice
- ma34 ACP -- Exercice
- Exercise: 3D point cloud
- Experience Data
- Calculations to find the characteristics of our cloud
- Exercise: 3D point cloud
- ma35 Système matriciel -- Exercices
- Vector Programming
- Partial Gaussian pivot method
- Choleski factorization
- Improve Jacobi
- ma40 Méthodes itératives
- Numerical simulation
- Iterative Methods
- Jacobi method
- Why does the 2nd case work?
- Calculation time
- ma41 Convergence de Jacobi avec inertie
- Add inertia to Jacobi
- Let's program inertia for Jacobi
- Let's study convergence
- Let's test other matrices with this algorithm
- Exercise 20.1
- Normalize
- Add inertia to Jacobi
- ma42 Surrelaxation pour Gauss-Seidel -- Exercice
- Exercise ma21
- Gauss-Seidel
- Gauss-Seidel overrelaxation
- Let's program overrelaxed Gauss-Seidel
- The good case
- Study by $w$
- Gauss-Seidel
- Exercise ma21
- ma50 Optimisation - Méthode du gradient
- Optimization problem
- Optimization problem with constraint
- The gradient method
- Study of the convergence of the gradient
- Optimization problem
- ma51 x.T A x sur un maillage en Numpy
- Let's calculate ${\bf x}^T \, A \, {\bf x} $ with Numpy
- Test case with A = 2 Id
- A real case
- Let's optimize
- Using a J function and a loop
- Using
np.tensordot
- Conclusion
- Let's calculate ${\bf x}^T \, A \, {\bf x} $ with Numpy
- ma52 Méthode du gradiant pour système matriciel
- Gradient et dérivée
- A x = b seen as an optimization problem
- Calculation of derivative
- Definition
- Calculate the derivative of J along a direction
- A symmetrical
- Gradient and derivative
- Calculation of derivative
- ma53 Notations du produit scalaire
- Writings of the scalar product
- ${\bf v} \,.\, {\bf w}$
- ${\bf v}^T \, {\bf w}$
- $<{\bf v}, {\bf w}>$
- Writings of the scalar product
- ma54 Gradient pour résoudre Ax = b -- Exercice
- The gradient method to solve A x = b
- Introduce inertia
- Optimal value of µ
- ma60 Méthode du gradient conjugué
- Conjugate gradient method
- Generate a base of $ℝ^n$
- The $A {\bf x} = {\bf b}$ case
- Calculation of $μ^k$
- Property
- 2nd attempt
- Let's work in the base of $\nabla J({\bf x}^i)$
- New calculation of μ
- Generate a base of $ℝ^n$
- Conjugate gradient method
- ma61 Système matriciel non linéaire
- Système matriciel non linéaire
- La méthode du point fixe
- La méthode du point fixe pour résoudre $A({\bf x}) \, {\bf x} = {\bf b}$
- Test
- Appliquons l'inertie
- La méthode de Newton-Raphson
- Système matriciel non linéaire
- ma62 Gradient conjugué -- Exercice
- Programmer le gradient conjugué
- Comparons avec le gradient simple
- Perfs
- Nombre d'iteration dans les 2 cas
- Un cas réel
- Comparaison gradient simple et conjugué
- Comparaison avec
lin.solve
de Scipy - Le gradient conjugué de Scipy (avec Lapack)